Broken symmetries and pattern formation in two-frequency forced Faraday waves 
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We exploit the presence of approximate (broken) symmetries to obtain general scaling laws gov- 
erning the process of pattern formation in weakly damped Faraday waves. Specifically, we consider 
a two-frequency forcing function and trace the effects of time translation, time reversal and Hamil- 
tonian structure for three illustrative examples: hexagons, two-mode superlattices, and two-mode 
rhomboids. By means of explicit parameter symmetries, we show how the size of various three-wave 
resonant interactions depends on the frequency ratio m:n and on the relative temporal phase of 
the two driving terms. These symmetry-based predictions are verified for numerically calculated 
coefficients, and help explain the results of recent experiments. 
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Symmetry arguments play a central role in our un- 
derstanding of pattern formation in nonlinear, nonequi- 
librium systems. They explain, for example, the ubiq- 
uity of simple patterns such as hexagons, squares and 
stripes observed in widely disparate systems, both in na- 
ture and in the laboratory Q. Group theory provides 
the natural language for describing and classifying pat- 
terns and, when combined with bifurcation theory, it de- 
termines which patterns may result from a symmetry- 
breaking instability It is by now standard practice 
to invoke symmetry arguments to determine which non- 
linear terms are present in the universal amplitude equa- 
tions that govern the pattern formation processes near 
onset. The detailed physics is manifest in these ampli- 
tude equations only through the numerical values of the 
coefficients of the nonlinear terms, yet the computation 
of these coefficients is often an arduous task that must 
be carried out numerically. 

In this letter we develop symmetry arguments in a 
somewhat nonstandard way, using parameter symmetries 
(symmetries which act on both amplitudes and param- 
eters; see e.g. ||) to determine the dependence of the 
most important coefficients in the amplitude equations 
on the parameters of the forcing function and on the 
damping. The specific problem we consider is that of 
gravity-capillary wave excitation on the free surface of a 
fluid subjected to the (vertical) forcing function 



F(t) = -(f m e im " t + f n e in " t + cc.) 

= \f m \ cos{muit + (j> m ) + cos(nojt + <p n ) 



(1) 



where m and n are relatively prime. This system, with its 
wealth of readily tuneable control parameters, has been 
a rich source of intriguing patterns j|, ||, [| [j], |8). Al- 
though a variety of nonlinear interactions influence the 
pattern selection process, resonant triads are especially 
important because they arise at second order in a weakly- 
nonlinear expansion and lead to a strong phase coupling 
between the three waves involved - they have been im- 



plicated in many Faraday wave pattern studies, both ex- 
perimental n d § 



10 . This letter 



and theoretical 
uses general symmetry arguments applied to several such 
three-wave interactions to shed some light on the role of 
the control parameters ( | f m | , |/„|, raw, nu>, 4> m , <j> n ) in 
the pattern formation process. 

Symmetry considerations, both spatial and temporal, 
have already gone a long way toward elucidating the gen- 
eral mechanisms by which different patterns arise. For 
example, when a single frequency lo is used, a sufficiently 
large forcing excites standing waves (SW) with primary 
frequency u/2. After one period T of the driving, the 
SW eigenmodes ipj are exactly out of phase and, due to 



the symmetry t — ► t + T, tpj 



-ijjj, are barred from 



quadratic (i.e. three-wave) interactions. The addition of 
a second commensurate frequency (when not a multiple 
of the first) breaks this particular discrete translation 
symmetry and replaces it with one on a longer timescale: 
t — » t + pT. If p is even this weaker symmetry constraint 
permits three-wave interactions by recasting the ipj as 
harmonic modes. Despite the success of this argument 
in explaining the emergence of hexagons when a second 
frequency is added Q it fails to predict either the size 
of the symmetry-breaking effect or its dependence on the 
control parameters. 

In this letter we take advantage of nearby (approxi- 
mate) symmetries to make further progress. The main 
idea is that, for weakly-damped small-amplitude dynam- 
ics, the real physical system (damped forced SW) is inti- 
mately related to an underlying system (undamped un- 
forced traveling waves (TW)) that is subject to more 
stringent symmetry requirements: invariance under ar- 
bitrary time translations t — ■+ t + t, and time reversal 
t — > —t. If the undamped but forced problem is Hamilto- 
nian the system is still further constrained. Although the 
addition of forcing and damping destroys these symme- 
tries, their influence persists. In fact, the broken symme- 
tries can be replaced by corresponding parameter symme- 
tries which remain intact for any damping and forcing. It 
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is only to obtain useful scaling laws that we must restrict 
ourselves to small damping and forcing. 

The overall strategy is illustrated in Fig. |l|. First, we 
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FIG. 1: Relationship of SW equations to underlying symme- 
tries. The origin denotes systems with time translation and 
time reversal symmetry. The undamped problem is Hamilto- 
nian. SW bifurcate from the flat state along the solid line, 
while the dashed vertical line illustrates a realistic experimen- 
tal path. 

express physical variables such as surface height h(x,t) 
in terms of the TW eigenfunctions of the undamped un- 
forced problem: 



N 



h(x,t) = J2J2 Z f^ e l{k '- x±n ^ 



(2) 



=1 ± 



In all cases N is finite (we consider only periodic pat- 
terns). Evolution equations for the TW amplitudes Zj 
must respect the appropriate spatial symmetries, which 
depend on the planform considered, and the following 
parameter symmetries: 



T T : Z 



± 
3 

± , 



"►-(*, 7), (fmJn) ^ (fmjn). (3) 



zf~zf, (t l7 ) 



Here 7 = 2vk 2 /uj is a dimensionless damping param- 
eter (y denotes viscosity and k(ui) is a characteristic 
wavenumber). Next, we assume the TW equations can 
be expanded in powers of 7, f m and /„ and so deter- 
mine their form for small damping and forcing. Lastly, 
these expanded TW equations are reduced at the pri- 
mary bifurcation through a center manifold reduction to 
yield the equations describing SW near onset. The proce- 
dure outlined above can be applied to numerous patterns 
and we shall look briefly at three: hexagons, two-mode 
superlattices Q], and two-mode rhomboids 0. We test 
the predicted dependence of SW coefficients by numeri- 
cally calculating them from the Zhang- Vihals model 
a quasipotential formulation of the Faraday problem ap- 
propriate for a deep layer of nearly inviscid fluid. This 
reduction (see |f0[ ) is done directly at the primary bifur- 
cation to SW and does not rely on the symmetry argu- 
ments we develop here. 

I. Simple hexagons. These common patterns are fa- 
vored by quadratic interactions requiring harmonic eigen- 



modes. We therefore take the coprime integer associ- 
ated with the instability (n, say) to be even, and set 
N = 3, Qi = fl 2 = ^3 = nuj/2 in the TW expan- 
sion (||). The neutral stability curves in the (k, f) plane 
( f 2 = \ f rn \ 2 + \fn\ 2 ) might resemble the sketch in Fig. ||b. 
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FIG. 2: Simplified sketch of when hexagons arise. With n 
even the modes excited by f n are harmonic (labeled "+") 
and those excited by f m are subharmonic ("—"). 

Spatial symmetries (translation, reflection, and 60° rota- 
tion) in combination with Eqs. (|^) lead to 

Z^ = vZ\ — p,f n Z^ + vd{Z\Zj2 + Z^Z^)Z^ 



0fm (Z 2 Z 3 



tifrnfn 
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^f,iz 2 z z , 



Z2Z3, 




+ z+Y,~4\zf\ 2 - 

The other five equations follow from symmetry. At each 
order in amplitude, up to third order, the leading order 
terms (in f m , /„, 7) have been kept. Terms of order 
\Z\ 3 f are neglected here (and in the examples that follow) 
because they produce no new qualitative effects. Note 
that the quadratic terms in Eqs. (0) split into two cases: 
m = I and m > 1. The time reversal symmetry k in (||) 
forces most of the coefficients to be imaginary at leading 
order. For example, we may write 

A = Mr 7 + iifii + Vm\fm\ 2 + Mn|/n| 2 + M 7 7 2 ) + ■•• ■ (5) 

Here /2 r ,/2i, etc. are real and, with an appropriate sign 
convention for the forcing term in the governing equa- 
tions, we can take jli > 0. The only coefficient which 
does not submit to a simple transcription of Eq. (|^) is v 
because, having chosen the wavevectors of the expansion 
to lie on the critical circle (itself a function of 7), we must 
have the detuning, Im(u), vanish at 7 = / = 0. 

At the bifurcation to SW, eigenvectors take the form 
Zf = e ±iv A 3 where tp = cj> n /2 - tt/4 + 0(7). A center 
manifold reduction yields 

A 1 = f iA 1 +aA 2 A 3 + A 1 (a\A 1 \ 2 +b\A 2 \ 2 +b\A 3 \ 2 ), (6) 

with equations for A 2 and A3 obtained by cyclic permu- 
tation. The normal form coefficients are given (to lowest 
order) by 



Mr, 



C 2 7 



(2f3i - 7i) sin$, m = 1 

&i\fn\~ ~ COS$, m>l 
\frn\ n ( C 3 COS 2 $, 



C17, 



C 4 ||/„| m - 3 sin 2 $, m>l /' (7) 
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Here $ = ir/4 — </>//2 with <f>f = mcf) n — n<f) m (this is 
the only linear combination of <j) m and <j) n which cannot 
be set arbitrarily just by redefining the origin in time). 
The Cj are known functions of v r , jl m , etc., indepen- 
dent of f m , /„, 7 (see p]). Note that a inherits the 
fmi fn dependence of the corresponding resonant term 
in Eqs. [| and, since at the bifurcation to SW we have 
|/n|j \fm\ ~ 7j it follows that a scales as 7™/ 2 if m = 1, 3 
and <-y(™+ m ~ 3 )/ 2 otherwise. Furthermore, a has a sim- 
ple harmonic dependence on the TV-invariant phase <E>. 
Both of these predictions are borne out in Fig. |3| for sev- 
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FIG. 3: Numerical calculations based on the Zhang Vifials 
equations (for a summary of the procedure see (a) |a(7)| 

for several choices of m:n jo. Solid lines show predicted 
scalings. (b) a(<£>) for m:n = 1:2 and 3:2, both with 7 ~ 0.001. 
These cases have identical scaling in 7 but are phase-shifted 
in $. (c) Range of stable hexagons in Eqs. (g|) for m:n = 5:4 
and (realistic) parameters such that 7 ~ 0.245. 



eral values of m:n. This figure also illustrates how the 
oscillations in a can lead to oscillations in the range of 
forcing amplitude over which hexagons are stable. Simi- 
lar oscillations have been seen experimentally j|, [| . One 
can understand such behavior as follows: while the lin- 
ear eigenmodes are phase-locked to the forcing /„ which 
parametrically excites them, the nonlinear forcing from 
the resonant terms (which involves / m ) attempts to in- 
duce a different temporal phase. These two influences 
may cooperate or compete, depending on the relative 
phase 4>f- 

II. Two- mode superlattices (2MS). These were 
observed with almost all forcing ratios by Arbell and 
Fineberg ^, ||] near the bi- critical point where the two 
modes driven by f m and f n onset simultaneously. The 
relevant three-wave interaction also includes a third lin- 
early damped mode (see Fig. |]b) with predominant oscil- 



(b) 





FIG. 4: Sketch of the linear problem for 2MS. We take m > n. 

lations at half the "difference frequency" . The 2MS states 
typically emerge from an underlying square or hexagonal 



lattice and involve multiple resonant triads |6|, g|. For 
simplicity, we focus here on a single characteristic triad 
(Fig. [|a), i.e., we treat an invariant subspace of the full 
2MS problem. A combination of spatial and temporal 
symmetries (^) again determines the form of the six TW 
evolution equations. In the vicinity of the bi-critical point 
these TW equations reduce to SW equations 



A 1 = XA 1 +A 1 (alA^ + blAil 2 ) 
A 2 = nA 2 + A 2 (clA,] 2 + d\A 2 \ 2 ) 



with coefficients of the form 
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(9) 



Here the Cj are known functions, independent of 7, f m , 
and /„ while &i, Pi, and the imaginary parts 

of the quadratic resonant coefficients in Z^, Z^, and 
, respectively. The most important effect of Hamilto- 
nian structure on the TW equations comes via these last 
three coefficients. Specifically, if at 7 = there exists 
a Hamiltonian Ti such that dZ^ 1 /dt = ^idTL / d Z^ , then 

<5j — Pi — Hi ■ The equality here is an overstatement since 
the equations (and hence the coefficients) can always be 
rescaled. One expects simply on ~ Pi ~ hi and thus, 
for small damping, 1 « d ~ — c. If, in addition, both 
primary bifurcations (to pure modes) within Eqs. (g) are 
supercritical then < —a ~ — d <C 1. Under these con- 
ditions stable 2MS states (mixed states within Eqs. (g)) 
inhabit the region shown in Fig. |^a. This region agrees 
quite well with the experiments of Arbell and Fineberg 
H despite being based on a simplified version of the 2MS 
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Driving Amplitude (g) 40 Hz. 



FIG. 5: (a) Region of stable 2MS under the assumptions of 
small damping, Hamiltonian structure, and supercriticality, 
(b) Experimental data, courtesy of J. Fineberg, reproduced 
from g. Since A ~ |/ n | - |/n| crit and fx ~ |/ m | - |/ m | crit the 
orientation of the axes in (a) and (b) is the same. 

states. The indifference of Eqs. (^) to m and n helps ex- 
plain why 2MS patterns were observed for virtually all 
forcing ratios 0. 

III. Two-mode rhomboids (2MR). These patterns 
are found H near the bi-critical point but, unlike 2MS, 
do not rely on damped modes; see Fig. ^. As before, 
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FIG. 6: Sketch of the linear picture for 2MR. We take m odd 
and n even. 



spatial symmetries and Eqs. (|J) determine the form of 
the six TW equations that reduce to SW equations at 
the bi-critical point: 

A x = XA 1 +aA 2 A 3 + A 1 (a\A 1 \ 2 +b\A 2 \ 2 +c\A 3 \ 2 ), 
A 3 = ^A 3 + (3A 1 A 2 +A 3 (d\A 1 \ 2 + d\A 2 \ 2 + e\A 3 \ 2 ), (10) 

with resonant coefficients of the form 



\.fn 



|/„| m 2 cos$. 



(11) 



Again <I> = 7r/4 — (pt/2 with 4>f = m<f> n — ruj) m . Hamil- 
tonian structure at 7 = implies on ~ Pi and thus 
diipi > 0, or equivalently, a/3 < (see |l3| for proof 
of an analogous relation in a class of self-adjoint hydro- 
dynamic problems). The condition a/3 < has enormous 
dynamical consequences for Eqs. ( J]~0| ) and is associated 
with the presence of modulated drifting waves and hete- 
roclinic cycles [jllj |l4|. We have confirmed the opposing 
signs of a and 0, as well as the predicted scaling \a\, 
1^1 ro /2 ^ f or a w [^ e ran g e f forcing ratios by 

numerically calculating these coefficients from the Zhang- 
Vihals equations. The exponential decrease of \a\ and 
|/3 1 with m + n helps explain why 2MR have been seen in 
jj], D for m:n = 3:2 and 5:4 but not with higher values of 
m + n: for m:n = 1:2 the same three-wave resonance is 
observed, but as part of a superlattice-II type state g. 

In conclusion, we have shown that with weak damp- 
ing the "broken" symmetries of time translation, time 
reversal, and Hamiltonian structure impose significant 
constraints on the equations describing SW near onset. 
These include scaling laws of resonant interactions with 
f m and /„, harmonic dependence on (f>j, and the dy- 
namically significant preference for resonant coefficients 
of opposite sign in the 2MR problem. We have verified 
numerically these predictions for the Zhang- Vihals equa- 
tions, and used them to interpret experimental results - 
they explain the oscillations seen with varying temporal 
phase in || and predict the correct region of param- 
eter space for 2MS states. We note, however, that it is 
not always clear whether the damping in a given experi- 
ment is "small enough" for our results to apply. This is 
particularly true when shallow containers are used and 
the damping due to boundary layer effects is important. 
It is unclear as well whether some procedures (such as 



the simple Taylor expansion in 7) would carry over to 
the full problem (Navier-Stokes equations with realistic 
boundary conditions) where the addition of damping con- 
stitutes a singular perturbation. 

Finally, we emphasize that the ideas developed in this 
letter can be fruitfully applied to a variety of other pat- 
terns, and could be extended to include more general 
periodic forcing functions. In particular, they indicate 
how one might control resonant wave interactions (and 
hence pattern formation) through a judicious choice of 
the harmonic content of the forcing function. Further- 
more, many aspects of our mathematical framework, de- 
veloped here for parametrically excited SW patterns, 
should carry over to other problems in which TW nat- 
urally occur, e.g. via Hopf bifurcation. It would be in- 
teresting, for instance, to apply our framework to the 
problem of SW pattern control via weak temporal mod- 
ulation of the control parameter Jl5|, . 
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